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I. INTRODUCTION 

Infrared instabilities near Fermi surface of degenerated 
Fermi gas are seen in many quantum systems. Supercon¬ 
ductivity/superfluidity is a well-known example, which 
appears not only in electron and atom systems but also 
in nuclear matter and quark matter [I|. Also the Kondo 
effect is caused by an infrared instability near the Fermi 
surface when interactions between degenerated fermions 
and sufficiently heavy impurities have non-Abelian na¬ 
ture (e.g., spin-flip interaction) d, d]. The Kondo ef¬ 
fect induces an enhancement in scattering amplitudes 
between a light fermion and a heavy impurity no mat¬ 
ter how weak an elementary coupling strength is. The 
Kondo effects are observed in experiments as a change 
of resistivity in a variety of systems from alloys to quan¬ 
tum dots. To investigate the Kondo effect, a number of 
theoretical methods has been developed: resummation 
of leading logarithm Q , the scaling-law approach Q , the 
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numerical renormalization group (RG) method Q, and 
so on. 

The Kondo effect occurs when a system is character¬ 
ized by the following four ingredients Q : (0) heavy im¬ 
purity, (1) Fermi surface, (2) quantum loop effect, and (3) 
non-Abelian interaction. These can be indeed identified 
in finite density QCD with dilute heavy quarks and thus 
the Kondo effect may take place, as recently suggested 
by one of the present authors Q. In Ref. emergence 
of the infrared instabilities was found in the scattering 
amplitude between a light fermion and a heavy impurity. 
The scatterings with impurities are governed by strong 
interaction involving several internal degrees of freedom 
responsible for the non-Abelian properties, such as an 
isospin for a D (B) meson in nuclear matter and the 
color for a charm (bottom) quark in quark matter. Note 
that, while a (pseudo-)spin flip plays a role in electron 
systems to induce a non-Abelian interaction, spin-flip in¬ 
teractions are suppressed in the large quark mass limit 
in QCD, leaving frozen spin degrees of freedom im. 
Nevertheless, QCD involves complex non-Abelian inter¬ 
actions associated with the isospin and the color degrees 
of freedom, and they give rise to intriguing phenomena 
in QCD matter. We call the Kondo effect whose non- 
Abelian property is provided by strong interaction the 
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“QCD Kondo effect.” When the Kondo effect exists in 
nuclear and quark systems, it will cause a large modifica¬ 
tion in the transport properties and excitations near the 
Fermi surface. Exploring the Kondo effect is thus useful 
to probe the nuclear and quark matter states. 

We investigate the quark matter composed of light 
(up, down and strange) quarks at high density with the 
heavy (charm or bottom) quarks embedded as impuri¬ 
ties. We show that a strength of interactions between 
the light-quark gas and heavy-quark impurity grows near 
the Fermi surface regardless of a weak coupling regime in 
high density QCD where the asymptotic freedom works. 
To this end, we deal with a free light-quark gas and 
dilute enough heavy-quark impurities, so mutual inter¬ 
actions between the heavy quarks are negligible. The 
non-Abelian interaction between the impurity and the 
scattering fermions, which is important in the Kondo ef¬ 
fect, is fulfilled by the color SU(n) symmetry (n = 3) 
causing the color-flip in the gluon-exchange interaction. 
The heavy and light quarks belong to the fundamental 
representation (3c), and a one-gluon exchange between 
two quarks i, j gives a product of color operators ti ■ tj 

with tk = {tl,... “^) being the generators of the color 

SU(n) group for the quark k. This operator gives the 
non-Abelian interaction between the heavy quark and the 
light quark. While the color SU(3) in quark systems is 
different from the conventional (pseudo-)spin SU(2) in 
electron systems, we will find that QCD governed by the 
color SU(3) nevertheless contains all the necessary con¬ 
ditions for the Kondo effect to emerge. 

It is important in the gluon-exchange interaction that 
the electric and magnetic components have different 
properties at finite density. The electric component gives 
the short-range interaction thanks to the screening ef¬ 
fect with the Debye mass in the long-range limit. The 
magnetic component, on the other hand, remains a long- 
range interaction, because it is affected only by the dy¬ 
namical screening M- In the present work, we include 
finite-range gluon exchanges with this difference in elec¬ 
tric and magnetic components taken into account. With 
this respect as well, we improve the preceding work by 
contact interactions Q. The unscreened magnetic com¬ 
ponent plays an important role in color superconductiv¬ 
ity, because it modifies a parametric dependence of the 
“energy gap” A on the QCD coupling constant g: the de¬ 
pendence changes from A oc (the screening case) 

to A oc e“'^ (the unscreening case) with some constants 
c and d [H, (see also Refs. Hill)- Thus, A is es¬ 
timated to be parametrically large, when effects of the 
magnetic component are taken into account. We also in¬ 
vestigate effects of the magnetic gluon exchange on the 
Kondo effect in the present work. 

In this work, we apply the RG method, which has been 
used for the color superconductivity [l^ [H, [I^ . We fol¬ 
low the perturbative calculation by assuming that the 
interaction coupling is small at the energy scale near the 
Fermi surface. Considering the scattering amplitude be¬ 
tween the heavy quark and the light quark up to one-loop 
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FIG. 1. Diagrammatic representation of a RG equation in 
a toy model. Thin and thick lines denote the propagators 
of light and heavy quarks, respectively. Grey blobs indicate 
scattering amplitudes at each energy scale. 

level, we will investigate how the scattering amplitude de¬ 
pends on the energy scale, and will estimate the Kondo 
scale in which the scattering amplitude becomes diver¬ 
gent and the system becomes a strongly coupled one. 

The article is organized as follows. In Sec. |lTl we 
present a toy-model analysis for the Kondo effect by us¬ 
ing the RG with a contact interaction. In Sec. IIIIl we 
describe the formalism of the gluon-exchange interaction 
in high density QCD and of heavy quark effective theory. 
We compute the scattering amplitude up to one-loop level 
and analyze the RG equation with the gluon exchange in¬ 
teraction. The final section is devoted to the discussions 
and the conclusion. 


II. KONDO EFFECT IN TOY MODEL 

For a simple illustration of the Kondo effect realized 
by color-exchange interactions, we briefly discuss a log¬ 
arithmic infrared divergence in the scattering amplitude 
at one-loop level 0 and then resummation of the lead¬ 
ing logarithm by a scaling approach Q. We compute a 
scattering process between an impurity (ff) and a scat¬ 
tering light fermion (-0) described by a simple contact- 
interaction Hamiltonian 

Hint = (1) 

c kl,ij 

where G (c = 1,..., — 1) are the generators of the in¬ 

ternal SU(n) symmetry Q. We keep a general integer n 
in the following analysis, but the symmetry group is sup¬ 
posed to be the color SU(3) group. Therefore, we will find 
that the Kondo effect is induced by the color-flipping in¬ 
teraction, although spin-flipping process, conventionally 
discussed in electroii systems, is suppressed in the large 
quark mass limit pll-ll3l|. The sign of G is chosen so that 
the interaction vertex at the tree-level gives an attrac¬ 
tion for the 3c channel and a repulsion for the 6c channel 
for n = 3. When the coupling is small, the perturbative 
calculation may be applied to the scattering amplitude. 
Both the impurity and the light fermion belong to the 
fundamental representation of the internal SU(n) sym¬ 
metry; k and I are the indices for the final and initial 
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states of the scattering light fermion, and i and j are 
those of the impurity {k,l,i,j = 1,..., n). 


A. Logarithmic infrared divergence 

To investigate an instability near the Fermi surface, we 
focus on scatterings of light fermions on the Fermi surface 
in the initial and final states. The scattering amplitude 
is shown in Fig. [T] up to one-loop level . At the tree level, 
the scattering amplitude is simply given by 


with a product of the color matrices 

Tkl,ij = )klit )ij- 


( 2 ) 


(3) 


At the one-loop order, the scattering amplitude between 
a light quark and a heavy-quark impurity 0 acquires a 
quantum loop effect as (see Appendix A for derivation) 




’ ^ J —E -I- le ki,ij j _ jg. 


P{E) 


(4) 


where the first and second terms come from the one- 
loop diagrams in the second row in Fig. [TJ respectively, 
and and noncommutative products of the 

color matrices on the interaction vertices 


’-{a) _ 




c,d k' i' 

= E . (6) 

c,d k' i' 

Infinitesimal complex displacements Eie specify bound¬ 
ary conditions as usual, and the first and second terms on 
the right-hand side correspond to the particle and hole 
propagations, respectively. E is the absolute value of the 
energy measured from the Fermi surface {E > 0). A den¬ 
sity of states p{E) at energy E can be expanded near the 
Fermi surface as p{E) ^ po + 0{E). The constant term 
po is finite in the presence of the Fermi surface, and thus 
the numerators of the both terms are finite in the infrared 
regime E ^ 0. Therefore, we find logarithmic infrared di¬ 
vergences as a Fermi surface effect. By decomposing the 
color matrices in Eqs. m and (O with the help of the 
identities. 


_ _ 

’^kl^ij o 


“ r, ( 4 „ r, ) Ekl, 


1 \ 1 
1 - I SkiSij - -Tki,ij, 


1 


1 




(7) 

( 8 ) 


^ We emphasize the four ingredients necessary for the realization 
of the Kondo effect by italic letters. 


we find that the scattering amplitude dU becomes 


(fc) 

ij 


„ n f dE 
= G^Po-Tki,,, I - 


dE 

E — ie 


E — ie 


(9) 


Here we find that, since the decomposed color matrices 
in Eqs. 0 and 0 contain the same terms, there is a 
large cancellation between the two terms in Eq. 0, i.e., 
the two one-loop diagrams in Eig. [TJ Only the last term 
in Eq. 0 survives to give the second line in Eq. 0. As 
mentioned above, this term has a logarithmic infrared 
divergence. If the matrices J2c(E)ki{E)ij hi the Hamil¬ 
tonian o were an Abelian type SkiSij, the infrared di¬ 
vergence in the one-loop calculation would completely 
cancel out and there were no infrared divergence. There¬ 
fore, the infrared divergence surviving in Eq. 0 reflects 
the non-Abelian properties of the generators F on the in¬ 
teraction vertices. This divergence appears with an arbi¬ 
trary SU(n) group, with the characteristics of the group 
seen in a simple factor of n in Eq. 0. 

In order to regularize the infrared divergence, we in¬ 
troduce an infrared cutoff parameter A. Then, we can 
rewrite the scattering amplitude 0 as 

= ( 10 ) 

Combining this with the tree-level amplitude 0 , we ob¬ 
tain the scattering amplitude as 



(5) 

+ 

o' 

II 

1 

)i'j ■ 

(6) 

— -I- G^ pQ--T};i ij / 

^ Je>a 


dE 

E — ie 


( 11 ) 


in the second term grows logarithmically in the infrared 
limit A —>■ 0, and becomes an important correction to the 
tree-level contribution in the first term. This work was 
extended to relativistic theories with the isospin/color 
symmetry in Ref. Q. While we focused on an illustra¬ 
tion of the Kondo effect by nonrelativistic expressions in 
this and the next subsections, we will proceed to the rel¬ 
ativistic theory on the basis of high-density QCD in the 
next section. 


B. Kondo scale from renormalization group 
equation 

Resummation of the leading logarithmic contributions 
can be handled by the renormalization group (RG) 
method. Following the description in Ref. [^, let us in¬ 
vestigate how the scattering amplitude Mkpij{A) behaves 
as the energy scale A changes. We find a RG equation 
for Mkpij{A) as 

Mki,ij (A — dA) 

dE , ^ 

( 12 ) 


Tt f 

= Mkpij{A) -\- G(A)^pQ—Tki,ij / 

^ Ja 


A-dA E — ie 
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Namely, the scattering amplitude in an energy scale 
A — dA is given by a sum of the tree-level contribution at 
A and the one-loop contribution integrated over a thin 
shell A — dA ^ A. Since this quantum correction has the 
same color-matrix structure as in the tree-level scatter¬ 
ing amplitude ([2|), one can absorb it into a renormalized 
coupling constant G'(A) giving a scattering amplitude in 
a lower energy scale 

(A - dA) = G(A - dA)Tfcz,y . (13) 


Following this renormalization, we obtain a flow equation 
for the effective coupling constant: 


d Ti 

A^G(A) = --G2(A)po, 


and its solution 

G(A) = 


G(Ao) 


l + ^G(Ao)log- 


(14) 


(15) 


Here Aq is an initial energy scale of the RG flow. We find 
that the running coupling constant G(A) has a Landau 
pole and diverges at A = Ak, where Ak is the so-called 
Kondo scale 


Ak = Ao exp (-) • (16) 

V npoG(Ao)/ 

The existence of the infrared divergence indicates that 
the system is a strongly interacting one near the Fermi 
surface, even though the interaction is weak in the ul¬ 
traviolet energy scale. The Kondo effect arises from four 
ingredients as we mentioned in the Introduction and ex¬ 
plicitly showed in this section. Namely, we found loga¬ 
rithmic terms originating from a loop integral in a scat¬ 
tering amplitude between a heavy-flavor impurity and a 
light quark carrying a Fermi momentum. Although these 
logarithmic terms cancel out if interactions are Abelian 
type, we showed that the cancellation is not perfect in 
the presence of non-Abelian interactions. Based on a RG 
equation, we found a Landau pole of the running coupling 
constant which becomes divergently large at the Kondo 
scale near the Fermi surface. 


III. QCD KONDO EFFECT 

In the last section, we have illustrated the emergence 
of the Kondo effect by using a contact-interaction model. 
We now investigate finite-range gluon interactions on the 
QCD basis, and how an effective four-quark vertex shown 
in Fig. [5] evolves as the relevant energy scale goes down 
to the infrared regime near the Fermi surface. 


A. Setup of problems 

We first look at a gluon propagator at high density. 
While the electric component is screened by the Debye 
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FIG. 2. An effective four-quark vertex induced by a scatter¬ 
ing between a light quark (thin line) and a heavy quark (thick 
line) through a gluon exchange interaction. 


mass mp, the magnetic component has only the dynami- 
cal screening effect called the Landau damping , lead¬ 
ing to unscreening in an infrared momentum region [H. 
An explicit form of each component is given by 


= < 


— Y for fi = V = 0 


m? 


i((5*J - q^q^) 


. 9o - l9p - '^ 1 '^ 


?,ko| 


(17) 

for /r, = i,j 


with g* = q'^/\q\, and the off-diagonal components 
are vanishing. Here we have taken the Landau gauge 
^ = 0 otherwise the gauge dependent term would exist 
^q^q'^/q'^. Since we will focus on a small momentum ex¬ 
change of the order of an infrared cutoff scale (A <C kp), 
the momentum dependence in the electric component 
can be ignored in the presence of the Debye mass. On 
the other hand, we need to maintain the momentum de¬ 
pendence in the magnetic component since the Landau 
damping term vanishes as the energy vanishes —>■ 0. 

We suppose that the mass of heavy quarks (charm and 
bottom) is much larger than the low energy scale in the 
system. A useful way to treat the heavy quark is to 
separate the hard scale momentum (of the order of the 
heavy-quark mass) and the residual momentum as 


= ruQV^ + k^, 


(18) 


with the heavy-quark mass toq and the four-velocity 
= (v°,v) = (y/l -I- |Fp, v) (gl-fl^. We note that 
the spatial component v in the four-velocity is different 
from the normal velocity u. Namely, v is related to u by 
V = u ■/x/l" — |up. The residual momentum k^ is sup¬ 
posed to be much smaller than the heavy-quark mass 
k^ <g; niQ so that the expansion with respect to l/mq 
is justified. A framework based on the 1 /mq expansion 
is called the heavy quark effective theory (HQET) [§- 
[T^ . which has been quite successful in describing heavy 
quark dynamics in a vacuum. While k^ could be at most 
of the order of AqcD in the vacuum, a heavy quark in a 
medium could typically receive a kick of the order of mo 
(a typical momentum transfer by one gluon exchange). 
Therefore, the validity of the HQET in a medium may 
be given by gp <C mq. 

The leading-order HQET Lagrangian reads 

-Chqet = Qviv-DQy + Aight + 0{l/mQ). (19) 









A positive-energy component of a quark field is denoted 
as Qy which was picked up by a projection 
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Q„(x)=e™«“-“i^Q(x), (20) 

where the hard momentum scale is factorized as a plane 
wave, and Qv{x) has only a residual momentum k^. 
While we examine a single heavy flavor below for sim¬ 
plicity, a generalization to multiple flavors is straightfor¬ 
ward. A covariant derivative + igA^ here is 

given by a derivative associated with the residual mo¬ 
mentum of the heavy quark a coupling constant g, 

and the gluon field with (c = 1,..., 8) 

being the generators of SU(n) group as in the toy model 
(El)- The QCD Lagrangian for light quarks and gluons 
is denoted as £iight- When the heavy quark has a finite 
three-dimensional velocity (|u | ^ 0) in the quark mat¬ 
ter, the magnetic gluons are coupled to the heavy quark 
dynamics at the leading order, because there is a non¬ 
vanishing coupling term between the spatial components 
V ■ Ain iv ■ D, while this coupling vanishes as the velocity 
vanishes (|iT| =0). 

As we already emphasized in the discussion on the toy 
model, the heavy-quark spin does not work as an ori¬ 
gin of non-Abelian interactions since the spin degrees of 
freedom are frozen in the heavy quark limit (toq —>■ oo). 
The Kondo effect is thus induced by the color-flip opera¬ 
tors This may be contrasted with the Kondo effect in 
condensed matter physics where the spin-flip interactions 
play a crucial role in the electron systems. 

We consider a light quark q scattering off a heavy- 
quark impurity Q at zero temperature, 

q\p)Q^ ^ q’^ip')Q\ (21) 

where a momentum exchange p' — p through the finite- 
range gluon interaction gives rise to a momentum depen¬ 
dence of the scattering amplitudes. Superscripts of the 
quark fields denote color indices. Let the scattering am¬ 
plitudes be G'^{p' —p) and G^{p' — p) as shown in Fig. [21 
and then effective four-quark vertices are represented by 

£eff = G^{p' - p)Y,q'^{p'h°v°{t<^rq^{p)Ql{ey^Qi 
a 

+gHp' -p)J 2 q'‘ip')Tvit‘^rq\p)Qlity^Ql 


( 22 ) 

We have divided the vertices into the temporal (T) and 
spatial (S) components, because the distinct properties of 
the electric and magnetic sectors in the gluon propagator 
(EZl) distinguish these components. By decomposing G"^ 
and G® into the partial waves, we define 

1 A 

GJ = - / d{cos 9) Pe{cos 8) G'^{p'- p), (23) 

Gf = - / d(cos0) P^(cos6l) G®(p'-p), (24) 

where 0 is a scattering angle between the initial and final 
states of the light quark. Since we focus on the S-wave 
scattering below, we suppress the subscripts of GJ^q and 
GLo as G'^ and G®, respectively. Then, the S-wave scat¬ 
tering amplitudes are expressed as 

-b G®?T-7) 0 ^^,(25) 

C 

where spinor indices of the gamma matrices are to be 
contracted with the light quark fields in the initial and 
final states. The heavy quark fields do not have the 
four-component spinor structure after the projections in 
Eq. (l2nil . and thus we do not have to take this structure 
into account above. 


B. Renormalization gronp eqnations 

With the setup in the previous subsection, we investi¬ 
gate the energy-scale dependence of the scattering ampli¬ 
tudes to derive RG equations for the effective four-quark 
vertices originating from the gluon-exchange interactions 
(fTTl) (see Fig. [2]). The scattering amplitudes are pertur- 
batively computed in the energy scale near the Fermi 
surface where the QCD coupling constant is sufficiently 
small at high density. We will show that the scattering 
amplitudes have infrared divergences, and that the RG 
flows evolve into regimes of strong couplings. 

The RG equation with gluon exchange interactions at 
one-loop order is depicted in Fig. |3l The tree-level scat¬ 
tering amplitude (the first line in Fig. |3|) reads 




= g^Tu,i 


dtPo[t) 


-1 

1 

m 




¥■1 


(po - - Ip - fc |2 - f|TO|)|po - fco|/|p - k I 


2 + 
D 


12 

23 


1 1 
'^3 


log 


8 -f i\ 




(26) 


where the kinematics is specified by 


bo - Pol = A, 


(27) 


Ip'I = Ip I = fcp, 

p' • p = cos 9 = t, 


(28) 

(29) 
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FIG. 3. Diagrammatic representation of a RG equation 
with gluon exchange interactions. Thin, thick and curly lines 
denote the propagators of light and heavy quarks, and gluons, 
respectively. 


where the second term is the Pauli blocking effect 
M- Note that the propagator of the heavy quark, 
i{l +‘^)/2v-{p — q), does not have the Pauli-blocking 
term, because the heavy quarks are treated as dilute im¬ 
purities. By applying the partial-wave decomposition to 
the gluon propagators D^i,[q — p) and D^v{q—k), we ob¬ 
tain the S-wave part of the scattering amplitude 
as 




I dqq^ 

9>fcp 


-Fi2(g)7° + -^2i(g)^ 


i + i! 


(34) 


where the integrand is given by 


with p' = p'/\p'\ and p = p/\p\- Here A is a small energy 
scale of the quantum fluctuation, as we are interested in 
an evolution of the RG equation near the Fermi surface. 
The three-dimensional momenta p' and p are approxi¬ 
mated to be the Fermi momentum fcp- In the second 
line in Eq. (1^ . we have defined A = (7r/2)mQA/A:p, 
which is a small quantity (A <C 1) when tod — nnd 
pL ~ k-p. We find that even the tree-level scattering am¬ 
plitude M% depends on the cutoff parameter A, be¬ 
cause of the energy-scale dependence A in the magnetic 
gluon propagator. The cutoff-parameter A is an energy 
measured from the Fermi surface. We assume that A is 
sufficiently smaller than kp, and hence that 0 ((A//cf)^) 
is neglected compared to 0{K/kp). We can identify tree 
level contributions to G'^ and G® from Eq. (1^ . Using 
expansion with respect to A ('Appendix IbI) . we find 


-^12(<Z) = (G^)' (35) 

+2vo\v\F2{q)G'^G^ + |F|2Fi(g) (GS)"| , 

and ^” 21 ( 9 ) with an interchange of the subscripts (1 O 2) 
in Eq. (1351) . The coefficient functions we compute below 
are 


Fi{q) = j dt 
F2{q)=J dt 


1 


voipo - qo) -v-p + \ v\\p\t + ie' 


t 


vo{po - qo) - V- p + \v\\p\t + ie 


(36) 


(37) 


In a similar way, we obtain the S-wave part of the scat¬ 
tering amplitude (the second one-loop diagram in 


9^ 



Eig.[3|) as 




(30) 




9 ' 2 1 

8 “h 

(31) 

^kl,ij ~ ^kl,ij 

[ dqq^ 

lq<k-p 

-^ 12 ( 9 ) 7 ° + ^ 21 ( 9 )^ 
kl. 

iX 





We again note that the unscreened magnetic gluon in¬ 
duces the A dependence of G® at the tree level, which 


was important in the color superconductivity [1^, [1 
Next, we proceed to computation of the quantum loop 


M 


( 1 ) 


= M 


(la) 


M 


( 16 ) 


effect at the one-loop level: 

where % and are the contributions from the 

first and second loop diagrams in Fig. [3] The scattering 
amplitude Ad ((“A is given by 




(a) 

kl,ij 


d^q 

Wy 

xiD, 


iS{q)^^v'^ 0 


*(1 -b it) 


2v{p-q) 

,{q - p)iDp„{q - k) . (32) 


We have the product of the color matrices as in 

Eq. dH), which will again play a role below as in the toy 
model calculation. The light-quark propagator at finite 
density and zero temperature is given by 


iS{q) = i 


i 

+ ie 


2Tr6{q'^)e{qo)9{kp - \q\) 


,(33) 


We note that the momentum integral in Eq. (l34)) is per¬ 
formed over the particle state above the Fermi surface 
due to the Pauli blocking effect in Eq. (1331) , while the mo¬ 
mentum integral in Eq. (I38|) is performed over the hole 
state below the Eermi surface. Since these integrals are 
divergent in the infrared energy region, we introduce in¬ 
frared cutoff parameters A > 0, and perform the integrals 
in Eqs. (|34)) and (l38l) up to g > A:f -b A and 9 < fcp — A, 
respectively. Eollowing from the integration, we obtain 
the energy-scale dependence of the scattering amplitude 
as 


=-^il2(A) + M£^(A). (39) 

Now that we have obtained the energy-scale depen¬ 
dence of the scattering amplitude, we can derive the RG 
equations as shown by using the toy model in Sec. |IT1 [see 
descriptions below Eq. mi Since the scattering ampli¬ 
tudes (1351) and ((351) have the spinor indices as well as the 
color indices, we compare the terms proportional to 7 *^ 






















7 


and V ■ 7 separately, and then obtain the RG equations 
for G'^(A) and G®(A), respectively as 

- = ^k\%ikF + AfTMkp + A) 

+-^J&(feF-A)2j-i2(/cF-A),(40) 
+ A)'J^2i(fcF + A) 
+^kihikF - A)2j-2i(fcF - A) 

where ^^ 12 ( 9 ) and J^ 2 i{q) are evaluated at q = kp ± A. 
The last term in Eq. m proportional to comes from 
the unscreened magnetic gluon exchange at the tree level. 
The last term in r.h.s. in Eq. dm was derived from the 
cutoff dependence in the tree level amplitude (I5T]) . Here 
we used the relation log(8+iA)/iA ~ log8/A+j(7r/2+n7r) 
for small A <C 1 with n being an integer, and left only 
the real part as the leading term in the infrared limit 
(A^O). 


C. Solution 


In order to obtain analytic solutions of the RG equa¬ 
tions (l40ll and (IdTI) , we shall examine the coefficient func¬ 
tions Fi{q) and ^ 2 ( 9 ) defined in Eqs. (1551) and (1571) . re¬ 
spectively. The kinematical variables are specified as 
Po = Ip I = ^F, I'c I = 'c, V’p = vk-p cos a, and go = I? I = 9- 
Then, we expand Fi{q) and F^iq) with respect to the 
nonrelativistic velocity of heavy impurities u <?; 1 , as 
shown in Appendix ICl Inserting Eqs. (IC5I) - (IC8I) obtained 
up to 0{v) into the RG equations (HUl) and (ITTl) . we find 
corresponding equations at 0(v) accuracy as 


dG^ 

dx 




2(2^)2 

+0{v^), 


n 
L X 


COS Qi - 1 / 'pv 2 

— —{J\f-2x)v (G^) 


(42) 


and 


dGS 

da; 




2(27r) 


—g'^g® - ^ 




(43) 


where a; = A/fcp is a normalized cutoff parameter, and a 
coefficient Af is defined by 


M = 


4 

G 




(44) 


with G being an eigenvalue of the Casimir operator 
)ki{F)ij. In case of the color SU(3) symmetry with 
n = 3, for example, the 3c channel gives G = —2/3, and 
the 6 c channel gives G = 1/3. When obtaining Eqs. (H71) 
and (H51) from Eqs. (HiH) and m, respectively, most of 


the terms cancel between the contributions from the two 
diagrams at the one-loop order as shown by the toy model 
[see Eq. ([9])] because of the common structures in the de¬ 
composed color matrices in Eqs. o and ([5]). 

By solving the coupled differential equations (1421) and 
(|43)) - '"^6 obtain a solution for G'’"(x) as 


G'^{x) = G'^(xo)^{x) 


+0{v^), 

where an initial value G"^(xo) of the 
scale xo = Ao/fep, and we introduced 

$(a;) = - 


n kp , X 


log-) 

cos a 

XoJ 



(45) 

is given 

at a 

5 

(46) 


1 -bGT(xo)- 77 ^ 3 ^ log — 


C{x) = 


2(27r) 


2 (27r)2 


Xq 


X So 


(47) 


We also obtain a solution for G®(s) as 
G®(a;) = -^Q - ^ (G^(so) -b Q) $(s) 

+ f { (1 + C(x))g^(xo) + 3G^(xo) - Q^‘^\x) 
+G(u), (48) 


where the terms originating from the Landau damping 
are proportional to 


Q 


2 g2 a; 

3 6*1'“^^;^ 


(49) 


Note that, since a factor of v was extracted in Eq. (1251) . 
G®(s) is G(l) when the scattering amplitude is 0{v) and 
is absent when the heavy impurities are at rest (s = 0 ). 

Interestingly, we find that <I> defined in Eq. (H5)) grows 
logarithmically as the infrared cutoff A approaches the 
Fermi surface. Therefore, both of the effective four-Fermi 
interactions G"^ (x) and G® (x) diverge at the Kondo scale 


/ 2(27r)2 \ 

which is smaller than an initial value Aq . When the scat¬ 
tering amplitude at the initial energy scale, say Aq ~ kp 
(xo — 1 ), is given by G"^(a;o = 1 ) — g^lk\, this diver¬ 
gence indeed emerges in the infrared regime since the 
Kondo scale is found to be much smaller than the Fermi 
momentum as 

( 2(27r)2\ 

Ak kp exp (- j <€. kp. (51) 


We note that the Kondo scale is common to several chan¬ 
nels of SU(n) symmetry, such as 3c and 6 c for n = 3. It 
is also straightforward to extend the above RG analysis 



















to a light-quark matter which contains heavy-antiquark 
impurities, and we find the same Kondo scale. 

The existence of the divergence in the infrared region 
indicates that the strengths of the four-quark vertices 
grow logarithmically, and the system goes into a strongly 
coupled one at low energy scale. The naive perturbation 
breaks down in such a regime as inferred by the emer¬ 
gence of the Kondo scale or the Landau pole. When 
n = 3, for example, the attractive and repulsive nature in 
3c and 6c channels, respectively, are strongly enhanced in 
the infrared regime. This result is not affected by the ex¬ 
istence of the long-range property in the magnetic gluon. 
It is also important to note that the infrared divergence 
appears not only in 0{v°) but also in 0{v^) in the scat¬ 
tering amplitudes and G® as shown in Eqs. (1^ and 

(H51) . 

Since effects of temperature smear the Fermi surface, 
we could regard the infrared cutoff energy near the Fermi 
surface as a temperature. In this sense, the physical in¬ 
terpretation of the Kondo scale Ak is regarded as the 
Kondo temperature at which the scattering amplitudes 
become divergent. This observation is consistent with 
the original analysis performed by Kondo in Ref. 0. 

IV. SUMMARY AND DISCUSSION 

In summary, we have shown that the Kondo effect is 
realized in quark matter with heavy quark impurities 
through the RG analysis of the scattering amplitudes be¬ 
tween a light quark and a heavy-flavor impurity. The 
emergence of the Kondo effect in quark matter is natural 
because this system possesses four necessary ingredients 
for the Kondo effects to occur, namely, heavy impurities, 
Fermi surface, quantum loop effect, and non-Abelian in¬ 
teractions. We have analyzed the scattering amplitudes 
with color-exchange interactions which are given by the 
(screened) electric and the (unscreened) magnetic gluon 
exchanges. By solving the RG equations for the effective 
interaction strengths and G®, we found that there ap¬ 
pears the Kondo scale at which the scattering amplitude 
diverges, and that the system becomes a strongly inter¬ 
acting one. We also found that the unscreened magnetic 
component in the gluon exchange does not contribute to 
the Kondo scale, but that it is important when a heavy 
impurity is moving at a nonzero velocity, which may be 
a more realistic situation in quark matter created in low- 
energy heavy-ion collisions. 

The Kondo effect is clearly a new idea in the field 
of QCD and will provide us with a novel viewpoint 
on the physics of high density quark matter. We 
conclude the paper by listing some of the interesting 
problems which we can investigate with this new picture. 

(I) Transport properties of quark matter at low temper¬ 
ature 

Increasing strength of the scattering amplitude im¬ 
plies a change of mobility of light quarks and thus 
will bring about significant change in the transport 


properties of quark matter. For example, similar to 
the ordinary Kondo effect in condensed matter, the 
QGD Kondo effect will increase the resistivity of light 
quarks in a quark-gluon plasma at low temperature. 
Also, we can expect that transport coefficients such as 
shear viscosity could be affected by the QGD Kondo 
effect. Furthermore, since the effect will be enhanced 
in proportion to the number density of charm and/or 
bottom quarks, we expect that anisotropic behavior in 
quark matter could be generated by inhomogeneous 
distribution of the heavy impurities. We will be able 
to study these problems experimentally in low-energy 
heavy-ion collisions at RHIG, GSTFAIR, and J-PARG 
which will create high density quark matter with charm 
quark impurities at moderate temperature. We also 
propose to study the Kondo effect and the modification 
of transport properties in lattice simulation for two-color 
QGD, which is free from the sign problem. Recall that 
the Kondo effect should occur for non-Abelian SU(n) 
groups including the two-color QCD with n = 2. 

(II) Kondo effect induced by heavy antiparticles 

It is straightforward to consider the Kondo effect when 
the heavy impurity is replaced by a heavy antiquark. 
We expect, however, that heavy quarks and heavy 
antiquarks play different roles in the Kondo effect 
because the matter is made of light quarks and the 
interaction between light quarks and heavy antiquarks 
should be different (we will discuss this point also in the 
next topic). In heavy-ion collisions, heavy quarks will 
be created together with heavy antiquarks. However, 
if heavy quarks and heavy antiquarks show different 
Kondo effects, we should be able to detect the effect 
of finite density or the existence of a quark Fermi surface. 

(III) Possible formation of the Kondo- Yosida state 
The strongly enhanced scattering amplitudes in the 
Kondo effect may suggest even the formation of a bound 
state of the light quark and the heavy quark. Such a 
bound state is discussed in the ordinary Kondo effect 
and is called the Kondo-Yosida (KY) state [l^. This is 
a bound state of the impurity and the light particle-hole 
pairs around the Fermi surface, which is induced by 
the enhanced interaction in an attractive channel (e.g., 
the spin-singlet channel in spin-spin interaction) in the 
Kondo effect. In quark matter, the heavy quark Q may 
lead to the formation of the color antitriplet (3c) KY 
state, given by the heavy quark and the light quark-hole 
pairs. In a similar way, the heavy antiquark Q may 
induce the formation of the color singlet (Ic) KY state. 
It may be interesting to note that the 3c KY state 
is colorful, while the Ic KY state is colorless. This 
difference may affect the hadronization properties of the 
heavy quark and antiquarks in a quark matter produced 
in the relativistic heavy ion collisions. 

(IV) Kondo effect in hadronic matter 

It is an important question to ask whether the Kondo 
effect can occur in hadronic matter (Table |T| . For the 
heavy quark Q, the heavy baryon Aq, which is the low- 
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TABLE I. The Kondo effect in quark matter and nuclear mat¬ 
ter for a heavy quark (Q) and a heavy antiquark (Q). The 
channels in which the Kondo effect occurs are checked by y/. 



quark matter 

nuclear matter 

heavy-quark 

Q V 

Aq 

heavy-antiquark 

Q 7 

Pq’ V 


est state qqQ baryon with spin-parity 1/2+ and isospin 
singlet {qq = ud), is the effective degree of freedom in 
nuclear mattefl Then, the non-Abelian interaction for 
Aq can be given by isospin- and spin-dependent inter¬ 
actions. However, there is indeed no such interaction, 
because Aq carries no isospin and the spin-flip interac¬ 
tion is suppressed by 1 /mg for large mass of heavy quark 
mg. Thus, we expect there will be no Kondo effect for 
the Aq baryon in nuclear matter. 

On the other hand, for the heavy antiquark Q, the 
heavy mesons Pg and Pq, where Pg (Pg) is the ground 

state of a heavy meson qQ with spin-parity 0“ (1~) 
and isospin doublet, exist as the effective degrees of 
freedom in nuclear matter 0. Interestingly, the interac¬ 
tion between a Pg (Pg) meson and a nucleon N has 
the isospin exchange induced by the isospin symmetry 
SU(2), and it can cause the Kondo effect. This was 
indeed demonstrated by using the contact interaction in 
the previous work 

(V) Possible continuity of Kondo- Yosida state from 
quark matter to hadronic matter 

In the previous two subsections, we have discussed that 
the heavy-quark (Q) can lead to the KY state in the 
quark matter, and not in the nuclear matter. On the 
other hand, the heavy-antiquark {Q) can lead to the KY 
state in both phases. We expect that the KY state for Q 
can have a continuous change, when the phase changes 
from the quark matter to the hadronic matter. This is 
analogous to the hadron-quark continuity discussed in 
the color-flavor locked (CFL) phase in color supercon¬ 
ductivity [I|. 

(VI) Competition between the Kondo effect and the 
color superconductivity 

It is interesting to consider the relationship between 
the Kondo effect and the color superconductivity. As 
analyzed in the present study, the Kondo effect is given 
by the existence of the ungapped Fermi surface. When 
the color superconductivity is realized, the Cooper pairs 
{qq) leads to the disappearance of the Fermi surface due 
to the energy gap, and hence the Kondo effect becomes 


^ We should recall that qQ meson in nuclear matter is an excited 
state, and is not considered in the ground state. 

® We note that and Pq are degenerate in mass in the heavy 
quark limit (utq oo), and hence both of them should be con¬ 
sidered. 


suppressed. Nevertheless, when the qQ interaction 
is much stronger than the qq interaction, the Kondo 
effects may overcome the color superconductivity. In the 
CFL phase, all the quark pairs {ud, ds, su) acquire the 
energy gap in three flavor case, and hence the Kondo 
effect will be competitive to CFL phase. In the two 
flavor case (2SC), however, only the ud pairs acquire 
the gap and ds and su pairs are not suffered from the 
Cooper pair formation, and hence the Kondo effect will 
be still realized for those free quarks. In any case, the 
competition between the Kondo effect and the color 
superconductivity will be an important subject in future 
study. 

(VII) Gauge invariance of the QCD Kondo effect 
All the calculation in Sec. nni was done in the Landau 
gauge ( = 0). Then there comes a natural question 
how the gauge dependence appears in the results, e.g., 
the Kondo scale. Let us first notice that the essence of 
the QCD Kondo effect can be captured by using the 
four-Fermi interaction as we showed in Sec. m implying 
that the gauge dependence in QCD calculation, if any, 
is expected to be a subdominant effect on the Kondo 
scale. In Sec. IIIII for the QCD calculation, one can 
explicitly show that at least the tree level contribution 
to the qQ scattering amplitude Eq. (E51) does not depend 
on the gauge parameter ^ (which is easily verified with 
the use of equations of motion for light quarks: the 
gauge dependent vertex factor u{p 2 )g'y'^f,q^tquu{pi) with 
9 + = P 2 ~ Pi vanishing due to {f) — m)u{p) = 0). 
The same argument does not hold for one-loop diagrams 
and thus it seems that there remains gauge-parameter 
dependence. Nevertheless we expect that the Kondo 
scale may be gauge invariant from the experience of the 
gap in color superconductivity. While the gap in color 
superconductivity apparently depends on the gauge 
parameter £ at higher orders, it is proven to be gauge 
invariant [2ll423l| . 

(VIII) Effects of self-energy 

In the present paper, we focused only on the effects in¬ 
duced by the presence of heavy quarks embedded in a 
light quark matter as impurities. However, it is pointed 
out that the quark matter itself becomes “non-Fermi liq¬ 
uid” due to logarithmic enhancement of the self-energy 
of light quarks within the framework of hard-dense loop 
approximation . Such a logarithmic enhancement 

may lead to a breakdown of perturbation theory at very 
low energies whose parametric dependence on the cou¬ 
pling constant is similar to the Kondo scale: exp(—a/g^). 
For the light-heavy quark scattering amplitude discussed 
in the present paper, the diagrams with insertions of 
the fermion’s self-energies appear at higher orders in the 
naive power counting. Although the origin of the loga¬ 
rithmic enhancement is different from that of the Kondo 
effect, one may need to include both effects for a com¬ 
plete description of the scattering amplitude, which will 
additionally require a wave-function renormalization for 
a light quark (logarithmic enhancement will be absent 
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for the heavy quark field included as a dilute impurity, 
since it does not form a Fermi surface and has a large 
mass scale regularizing the infrared dynamics). It is in¬ 
teresting enough to investigate how the Kondo effect is 
modified by the inclusion of non-Fermi liquid effects of a 
light quark matter. 
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Appendix A: One loop contribution in toy model 

We derive Eq. (|3]) in Sec.|lT]in the nonrelativistic limit 
(see Ref. Q and Sec. uni for relativistic computations). 
When the heavy impurity has a finite mass M and a 
momentum P, the loop integral is given by (see Fig. a 




2 cjria) 


d^k 


1 




1 1 1^12 1 1 ^, 2 ^ 1 1 ^^-. r |2 

iTrb r + r - r + + P - * r' 


2m 


2M' 


2m 


2M" 


te 




d^k 


1 




(Al) 


where the initial and final light quarks with mass m have 
the three-dimensional momenta p and p' on the Fermi 
surface (|p| = \p'\ = /cp), and k is the three-dimensional 
momentum in the loop integrals (cf. Ref. Q). The first 
integral on the right-hand side gives the contribution 
from quarks above the Fermi surface, and the second one 
gives that from holes below the Fermi surface. When the 
heavy impurity mass M is sufficiently large, the terms 
proportional to 1/M are neglected, and thus Eq. (lAll) is 
reduced to 


M, 


( 1 ) 

kl,ij 




(a) 


kltij 


|fc |>fcF 


d^k 

(^ 


ep - —|fc p -I- ie 
2m 


(A2) 




’\k |<feF 


d^k 


CF - ^\k p - ie 


with the Fermi energy ep = fc|/2m. By using the abso¬ 
lute value of the energy measured from the Fermi surface, 




for particles (|ft | > fcp), and 


(A3) 




(A4) 


for holes (|fc | < kp), Eq. (IA2I) is rewritten as 

“ PiE) 






0 —E is 

p{E) 


E — ie 


dE 

dA, 


(A5) 



FIG. 4. One loop diagrams in a toy model. 


with a density of states p{E) = i^rrif‘G\JE E ep/(47r2). 
As a result, we obtain Eq. Q. We can approximate p(E) 
as a constant po = (2m)^/^y^/(47r^) for small E <C ep. 


Appendix B: Integral formulas 


The integral appearing in the loop calculations can be 
carried out as 




where 


h(t) 


1 

2(1 - t)-k A/v'2(l - i) ' 


(Bl) 


(B2) 


To examine the infrared singularity arising near the 
Fermi surface, we may expand the following integrals 
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with respect to |A| <C 1 as 


After some calculation, we obtain Fi (q) as 


J t'^h{t)dt = -l + ^\ogj + --- , (B3) 

J ^{l-f)h{t)dt = l-+ ■■■ , (B4) 

1 

+ 5 + 

/■^ 2 .2/ X , 27r ,_2/c, 17 1,8 

JJ = +24-3“>®I + ’-' 

(B6) 

J ^{1 - f) h^{t) dt = -^ + ^logj + ■■ ■ . (B7) 


Appendix C: Expansion by a small velocity 


Substituting the kinematics specified in the beginning 
of Sec. IIII Cl the coefficient functions Fi{q) and F 2 {q) in 
Eps. (131)11 and dSZl) are rewritten by 


^^ 1 ( 9 ) 

F2{q) 


1 


I vo{kF — q) — kpv cos a + vqt + ie ’ 


dt—-- 7 - - -. 

I vo[kF — q) — kpv cos a + vqt + le 


(Cl) 

(C2) 


Fi{q) = 


vq \ 


log 


VT+lA — V 

y/Td-v^ + V 


(C3) 


-log 


q — (Vl + + r;)(\/l + v'^ — v cos a)kF\ 

q — (y/Td-^ — v){y/Td-^ — v COS a)kF I ’ 


and ^ 2 ( 9 ) as 


F 2 {q) = —^2 - {voikp - q) - -cgcosa) Fi(g)|, (C4) 

when |vTT ^((7 — fcp) +'yfcr cosa| > vq, because the 
energy scale Aq is the starting energy scale in the RG 
equation. We find series representation of Fi{q) and 
7^2 (g) at a small velocity u <C 1 as 

fi(g = fcF+A) = --+ v + Oiv^), (C5) 

F2iq = kF+A) = ^^^^ + Oiv^), (C6) 

when the light-quark momentum q is larger than the 
Fermi momentum, and 

, , 2 2fcFCOsa , , 

fi(g = fcF-A) = -+ v + 0{v^), (C7) 

F2iq = kF-A) = ^^^v + 0{v% (C8) 

when the light-quark momentum q is smaller than the 
Fermi momentum. 
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